Dynamics of grain ejection by sphere impact on a granular bed 



00 

o 
o 

(N 
o 

Q 



S. Dcbocufi^^ P. Gondrct\ and M. Rabaudi 
^Umv. Pans-Sud, Umv. Pans 6, CNRS, Lab. FAST, 
UMR 7608, Bat. 502, Campus Univ, 91405 Orsay, France and 
''Univ. Paris 6, Univ. Paris 7, CNRS, Lab. de Physique Statistique de I'Ecole Normale Superieure, 
UMR 8550, 24 rue Lhomond, 75231 Paris Cedex 05, France 
(Dated: December 17, 2008) 

The dynamics of grain ejection consecutive to a sphere impacting a granular material is inves- 
tigated experimentally and the variations of the characteristics of grain ejection with the control 
parameters are quantitatively studied. The time evolution of the corona formed by the ejected 
grains is reported, mainly in terms of its diameter and height, and favourably compared with a 
simple ballistic model. A key characteristic of the granular corona is that the angle formed by its 
edge with the horizontal granular surface remains constant during the ejection process, which again 
can be reproduced by the ballistic model. The number and the kinetic energy of the ejected grains 
is evaluated and allows for the calculation of an effective restitution coefficient characterizing the 
complex collision process between the impacting sphere and the fine granular target. The effective 
restitution coefficient is found to be constant when varying the control parameters. 

PACS numbers: 45.70.-n,45.50.-j,83.80.Fg,96.15.Qr 
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I. INTRODUCTION 



Impact cratering has been recognized as an important ge- 
ologic process for the last decades when the lunar craters 
have been finally attributed to impact structures rather 
than giant volcanoes as believed until 1950's The 
planetary impact craters such as the ones observed com- 
monly on the Moon or the Earth result from very high 
energy impacts of meteorites and thus involve numerous 
and very complex phenomena such as shock and rarefac- 
tion wave propagation, melt and vaporization of the pro- 
jectile and target materials, together with excavation by 
displacement and ejection of the target material [l|. In 
light of evidences of the discrete and sandy nature of 
planet surface, laboratory scale experiments of high en- 
ergy impacts on granular materials were conducted 0, Q ■ 
Since a few years, physicists have conducted laboratory 
scale experiments with rather low impact energy on gran- 
ular matter, interesting in the crater morphology and 
searching for scaling laws for the crater size [1, H S 13 ■ 
Even though their energies are typically many orders of 
magnitude smaller than those of meteorite impacts, these 
small scale experiments on granular impacts may be rel- 
evant to planetary impact processes, as the progression 
of crater morphologies as a function of impact energy 
has been shown to mirror that seen in lunar craters [j]. 
In these impact experiments, physicists have also been 
interested in the penetration of the imp acting sphere in 
the granular target [1, d, O, [O, [H, El, El ■ Indeed, de- 
spite recent progress on the complex rheology of granular 
matter [l5| . the penetration dynamics of a solid sphere 
into a granular medium is still difficult to understand 
well as it involves both the complex drag resulting from 
frictional and coUisional processes, and the final stop in- 
volving the complex "liquid/solid" transition exhibited 
by granular matter [l^. The penetration dynamics of 
the impacting sphere and the grain ejection have been 



shown to be very different when the granular material is 
not dense but loose: A spectacular thin granular jet can 
raise very high after the impact as first demonstrated 
in Refs [2, E3j EB|- The effects on this granular jet of 
the interstitial fiuid [l^l and of the initial packing 
fraction of the target [21| have then been studied. In the 
dense case, no granular jet but a growing granular corona 
is seen after the impact. These different kinds of grain 
ejection can be related to similar kinds of liquid ejec- 
tion consecutive to the impact of a droplet into a deep 
or thin layer of liquid that have been first filmed by j22l | 
and then studied extensively [H, [23|. Much less stud- 
ies focusing on the grain ejection have been performed 
in the dense granular case: Ogale et al. [2^ have mea- 
sured the mass of the spilled-over grains and Boudet et 
al. [26| have proposed a model of ejection from the crater 
growth in a layer of thickness small compared to the size 
of both projectile and target grains at low impact veloc- 
ities {r^ 1 m/s), while Cintala et al. [l^l measured ejec- 
tion speeds and angles of grains at high impact velocities 
(~ 1 km/s). Another type of experiments concern the 
impact of one bead with a granular target made of the 
same beads (28j . Such impacts and grain ejections by 
an impacting projectile have been recently simulated in 
limiting cases, for which the impact energy is very low 
or very high [29|, l30l| or for which the projectile size is 
about the grain size jsil, [13, HI] . We focus in the present 
paper on the dynamics of the granular corona formed by 
the ejected grains from a dense and deep pile upon low 
speed (~ 1 m/s) impacts. In section |lTl the experimen- 
tal setup is described together with the measurements. 
A simple ballistic model is then presented in section IIIII 
and compared to the experimental data. The results are 
discussed in section IIVI and section |V] ends this paper 
with a conclusion. 
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II. EXPERIMENTS 



Each experimental run consists of dropping a solid sphere 
into a granular medium. Four different steel spheres of 
density ps=7800 kg/m^ are used as impactors, with dif- 
ferent radius R = 5.15,6.75,7.55,9.50 mm and masses 
M ranging thus from 4.5 g to 30 g. The steel sphere is 
initially hold by a magnet through a semi-spherical hole, 
so that the sphere can be dropped without any transla- 
tional nor rotational velocity by pulling up the magnet. 
The sphere is released directly above the center of a con- 
tainer and falls along its axis. Being dropped from the 
height h above the granular surface which is varied be- 
tween 8 cm and 60 cm, the sphere thus impacts the gran- 
ular material with the velocity Uc varying from 1 m/s 
to 4 m/s and energy Ec = Mgh ranging from 3.10^'^ J 
to 2.10~^ J. The target material consists in sieved glass 
beads of density pg = 2500 kg/m^ {Ps/Pg = 3.1) and 
mean diameter 2r = 0.4 ± 0.1 mm {R/r ^ 25 — 50), thus 
of mass m = 8.10~-' g {M/m - lO'' - 10^). The granular 
material fills the cylindrical container of diameter 19 cm 
and height 26 cm. The size ratio of the container diam- 
eter over the sphere diameter is always larger than 10, 
so that there is no influence of the radial confinement by 
the lateral walls of the container neither of the bottom 
wall [l^ as the height of the packing is large. Before each 
drop, the granular medium is prepared by gently stirring 
the grains with a thin rod. The container is then over- 
filled and the surface levelled using a straightedge. The 
typical value of the solid volume fraction of the packing is 
60%. Each impact experiment is lighted from the front 
to enhance contrast between the grains and the black 
background, and side view images are recorded by a high 
speed camera at the rate of 500 imagcs/s and resolution 
of 0.16 mm/pixel. The 256 gray level images are thresh- 
olded to identify grains. 

Figure [1] shows a sequence of side view images separated 
by 30 ms illustrating the ejection dynamics of the grains 
after the sphere impact defined as the time t = 0. One 
can see few isolated ejected grains above a dense corona of 
grains that expands radially and vertically. The amount 
of ejected grains can be quantified by measuring the ap- 
parent surface area Atot of all the grains in each image. 
The evolution of Atot as a function of time is reported in 
Fig. [5] for several experiments corresponding to different 
dropping heights h of the same sphere. After the impact 
at t = defined as the first contact between the bottom 
of the sphere and the granular surface, Atot increases up 
to a maximal value denoted Atot max at the time 
before decreasing. Each curve corresponds to a single 
impact experiment, thus without any ensemble averag- 
ing, but with a little smoothing by a slide-average over 
a time window of 10 ms so that no data appears before 
t = 4 ms in the reported figures. Increasing the drop- 
ping height h of a given sphere, the values Atot max and 
tAtatma:! iucreasc, accounting for the increase of both the 
number of ejected grains and the dynamics duration. The 
same goes when keeping constant the dropping height h 




FIG. 1: Sequence of side view images separated by 30 ms, 
sliowing the typical time evolution of ejected grains after the 
impact of a steel sphere of radius 7? = 6.75 mm dropped from 
the height /i = 30 cm {Uc = 2.45 m/s, = 30 mj). 

and increasing the mass M of the impacting sphere. It 
is worthnoting that the exact relation between the mea- 
sured apparent surface area Atot of the ejected grains and 
the real total number of ejected grains N^^j is not straight- 
forward. For dilute zones of ejected grains, all the grains 
can be seen in the 2D images, but there is still the prob- 
lem of focusing: A far grain appears smaller than a close 
grain so that the grains do not have the same apparent 
area. The most problematic case concerns however the 
dense zones where some grains can be hidden by other 
grains. Using the fact that the ejection process is ax- 
isymmctric, the number of ejected grains N^j is related 
to the measured apparent surface area Atot by: 

Nej =nW'K Atot max, (1) 

where w is the unknown corona thickness and n the un- 
known number of grains per unit volume. The unknown 
product nw will be deduced in the following by compar- 
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FIG. 2: (Color online) Time evolution of the total apparent 
surface of ejected grains Atot for an impacting steel sphere of 
radius R — 7.55 mm dropped from different heights h = 13, 
23, 33, 43, 48 and 58 cm (from bottom to top). 




FIG. 3: (Color online) Geometrical parameters characterizing 
the corona of ejected grains. 



ing our scaling laws for the grain ejection with scaling 
laws already known in the literature for the crater size. 
In impact experiments of a drop onto a liquid layer, a 
beautiful corona is classically observed [H, US] . In such 
a liquid case, the corona is easy to define and extract as 
the liquid is a continuous medium. In the granular case, 
the corona is less easy to define as the ejected grains are 
individual entities. The granular corona is here defined as 
the largest connected part of grains in each side view im- 
age. We have checked that the size of the corona does not 
depend significantly on the lighting and contrast of the 
images. To investigate the time evolution of the corona 
and characterize its shape which is basically axisymmet- 
rical, its contour is extracted for each image. An example 
of such a contour is drawn in Fig. [31 Note that this cor- 
responds to the external contour of the ejecta curtain. 
This allows to define the bottom diameter Df, and the 
top diameter Dt as the minimal and maximal diameter 
of the corona. The height H of the corona is measured 
as the distance between the mean vertical position of the 
corona top contour and the granular surface level before 
impact. As the corona lateral edge appears quite straight 
except in a small zone at the base of the corona, we ex- 
tract also the angle Q formed by the corona edge with the 
horizontal, by a linear fit of the straight portion. 
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FIG. 4: (Color online) Time evolution of a) the height H, b) 
the top diameter Dt, c) the bottom diameter D;,, and d) the 
angle of the corona Q for the same experimental parameters as 
in Fig. [2] The horizontal line in d) is the time and ensemble 
average of Q(t) for the 65 experiments. 



4 



The time evolution of the corona, in terms of its height 
H, its top and bottom diameters Dt and I?f,, and its edge 
slope 0, is displayed in Fig.|4]for the same experiments as 
in Fig. O The expansion of the corona is demonstrated 
in Fig. [4^ by the increase of its height up to a maximal 
value denoted Hmax at time i_f/„„^ before decreasing. For 
a given sphere, Hmax and time tH^^a: increase monoti- 
cally with the dropping height h, so that the different 
curves of Fig. [4^ appear in order. Note that for large im- 
pact energies (large impact heights h), H decreases to a 
significant non-zero final value, because of the final crater 
rims lying above the initial free surface In the same 
time, the top diameter Dt and the bottom diameter Db 
increase with time, as shown in Fig.|4)3 and c, up to their 
maximal values, Dtmax and Dhmax, when the corona dis- 
appears and its height vanishes. The evolutions of Dt and 
Di, arc different: About linear for Dt but parabolic for 
Dh. Besides, Dt increases significantly with the dropping 
height h whereas Df, does not vary so much. Note that at 
the nearly end of the corona life, for vanishing height H, 
the values of Dt and Db become very noisy and are thus 
not shown here. The angle 6 the corona edge forms with 
the horizontal, shown in Fig. [3Jl, is roughly constant as 
a function of time and whatever the dropping height h, 
and equal to about 56° with relative variations of ~ 5%. 
The same kind of evolutions of all these parameters are 
observed for the different tested spheres. 

The corona evolution reported here for granular impacts 
can be compared to the corona evolution for the liquid 
case [13]. For the liquid case, the corona base spreads 
radially as the square root of time whereas it seems more 
close to a quadratic evolution for the granular case (see 
Dii{t) in Fig. Hj:). The scaling laws for spreading is thus 
different in the granular and liquid cases, with inverse 
curvatures in the Di,{t) plot. Besides, the angle of the 
granular corona is found independent of the impact ve- 
locities, as found for the liquid corona, with a value closer 
to the case of shallow liquid layers than the case of deep 
liquid layers [1^. This may be related to the coUisional 
chains that redirected the impact velocity within a few 
grains layer only [33 |. 

With the measurements of the geometrical parameters 
of the corona, we can now define the apparent area of 
the corona Acor, that is related to the amount of ejected 
grains contained in the corona, and compare it to the 
total apparent area of ejected grains Atot, related to the 
total amount of the ejected grains, both in the corona 
and isolated. Acor is measured as the area included in- 
side the corona contour, which is not far from the area 
H{Dt -t- Db)l2 corresponding to the approximate corona 
trapezium shape in the images. The time evolution of 
Acor is reported in Fig. [S] for the same experimental 
parameters as in Figs. 2 and 4. Even if A^or is al- 
ways smaller than Atoti the evolution of Acor is quali- 
tatively the same as the one of Atot (Fig- [l]). This is 
confirmed by the quantitative comparison of the corre- 
sponding coordinates of the curve maxima, (tAt<,t„„^, 




FIG. 5: (Color online) Time evolution of the apparent surface 
area of the corona Acor for the same experimental parameters 
as in Fig. [2] 
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FIG. 6: (Color online) a) Time fAcor,nax for maximal value of 
the corona apparent area A^or as a function of time tAtotmo.^ 
for maximal value of the total apparent area Atot- The solid 
line is of slope 1. b) Maximal values Acorn%a:c of Acor as a 
function of maximal values Atot max of Atot- The solid line 
is a linear fit. Data symbols are for the 65 experiments with 
difi'erent dropping heights h for impacting spheres of radius 
R = 5.15 (o), 6.75 (□), 7.55 (o) and 9.50 (V) mm. 



Fig. [S] for all the experiments (all dropping heights, all 
different impacting spheres). The evolutions of ^cor and 
Atot are synchronized in time as illustrated by the equal- 
ity *A,<,,„„, ^ *At„t„„x (Fig. [6li). Besides, A^ormax is 
found proportional to Atot max with the same ratio for 
all experiments: A^ormax — 0-75Atotmax (Fig. db)- AH 
this suggests that the investigation of the corona dynam- 
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ics is a good first order for the study of the dynamics of 
grain ejection due to an impact. Let us now interpret 
these experimental results by a simple ballistic model. 



III. BALLISTIC MODEL FOR GRAIN 
EJECTION 

The dynamics of the grains is found axisymmetric 
and rapidly contact forces between grains play no role 
(Fig. [1]). Furthermore air friction can be estimated 
and is negligible compared to the grain weight. Thus 
rapidly after impact, each grain trajectory corresponds to 
a parabolic free flight under the action of gravity alone. 
Using these preliminary remarks, we now build a very 
simple axisymmetric model that reproduces the observed 
corona dynamics in the (r, z) plane of cylindrical coordi- 
nates. 

As a first attempt, we assume that all grains start at 
the same time t = 0, from the same position (ro,0), in 
the same initial direction making an angle a with the 
horizontal, but with different velocity amplitudes. It is 
thus easy to show (Fig. [7^) that at any time t all the 
grains will be located on a cone making the same angle 
a with a bottom diameter 2ri increasing quadratically 
with time as 2ri(t) = 2ro + gt^ / tana. Furthermore, if 
Mo is the maximal amplitude of initial velocity, the top 
diameter of the cone 2r2 increases linearly with time as 
2r2 (t) = 2ro + 2uo cos a t and the height of the cone Z2 
evolves quadratically in time as Z2(t) = uq sinat — gt^ /2. 
These analytic results are quite similar to the dynam- 
ics observed in Fig. [1] and to the time evolution of the 
corona parameters H, Dt, Dj, and 6 presented in Fig. H] 
The hypothesis of an instantaneous release of all the 
grains at the same place and the same time t = is 
however clearly oversimplified. Indeed grains are ejected 
during a time interval St of the order of the penetra- 
tion time R/Uc ^ 5 ms and, as the projectile decelerates 
the ejection velocity of the grains should be a decreas- 
ing function of time. In Fig. [7)3, we plot the apparent 
cone angle 0{t) of the moving grains when released with 
the same ejection angle a but with various models of 
time decreasing velocities uo{t) during the time interval 
6t = 5 ms: A linear decrease [uej{t) = uo{l — t/5t)], a 
quadratic decrease [uej{t) = uo{l — t^/St'^)], or an expo- 
nential decrease [uej{t) — uq eicp{—t / St)] from the initial 
value Mo = 1 m/s. Except at very short times the edge 
angle 9 remains constant and always close to the ejection 
angle a. While the time evolutions for the top diameter 
2r2 and the height Z2 of the apparent cone are obviously 
the same than for an instantaneous release, it appears 
to be also the same but with a small delay time St for 
the bottom diameter 2ri. As this delayed grain ejec- 
tion and initial velocity decrease have no visible effect 
in the corona dynamics, we will use in the following the 
first simpler model where all the grains are ejected at the 
same time with ejection velocities < Uej ^ uq. 
In the experiments we measured the time evolution of 
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FIG. 7: (Color online) a) Ballistic model for grains ejected in- 
stantaneously from the same position ro in the same direction 
making the angle a with the horizontal, with different veloc- 
ity amplitudes. Trajectories are drawn in dashed lines. At 
any time grains align along a straight line forming the corona 
edge of angle = a drawn in continuous lines, b) Time evo- 
lution of the angle of the cone edge if grains are ejected 
during a finite time 5t = 5 ms in the direction a = 53°, with 
a linear, quadratic or exponential decrease of ejection velocity 
from the initial value uo = 1 m/s. 



the corona parameters, i.e. the optically opaque zone 
formed by ejected grains, when viewed from the side. 
The relation between this corona that corresponds to a 
2D projection of the 3D real dynamics with this latter 
is not straightforward. But the bottom diameter Df, and 
the edge angle 6 are clearly the same. Furthermore, Fig. [HI 
shows that most of the grains are inside the corona. In 
the following we will then assume that the height H and 
the top diameter Dt of the corona in the experiments are 
identical to the height Z2 and the top diameter 2r2 of 
the cone in the ballistic model. Assuming this, the time 
evolution of H, Dt, Di, and 6 arc given by the following 
set of equations: 

H{t) = uosinat - gt'^/2 + H{0), (2a) 

Dt{t) = DtiO) + 2uq cos at, (2b) 

Db{t) = DbiO)+gt'^/tana, (2c) 

6 ^ a, (2d) 

An initial height H{0) has been introduced in the dy- 
namics of H(t) as we observe in the experiment an initial 
swelling of the granular bed before any grain ejection. 
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This swelling corresponds to a deformation of the sub- 
strate when the projectile starts to penetrate the bed. 
This deformation process is clearly not contained in our 
ballistic description. The horizontal position of ejection 
To is related to the projectile size: We checked that Dt{0) 
and Db{0) increase with the impacting sphere diameter 
2R. The set of equations can be written in 

dimensionlcss form as: 



H„ 



H{t)-H{0) 
-H{Q) 

Dtjt) 

A(0) 
Db{t) 
DM 



= 2 



t 



t 



2uo cos a 



= 1 



A(0) 
9 



Db{0) tana 



, (3a) 
(3b) 
(3c) 

Mo sin a / g and 



with the two scaling parameters: 

For each of the 65 experiments, our experimental data 
for the height H{t) and the top diameter Dt{t) of the 
corona are well fitted by the ballistic model equations 
as shown in the non-dimensional plots of Fig. [5^ and b. 
The evolution of H is basically quadratic, while Dt is 
rather linear in time. The fitting parameters mq sin a for 
H{t) and uocosa for Dt give access to the typical grain 
ejection velocity ug and to the typical ejection angle a 
for each impact experiment. The time evolution of the 
bottom diameter of the corona Df){t) is shown in the non- 
dimensional plot of Fig. [Ht, with the already deduced 
fitting parameters: The predicted quadratic increase of 
Dh is rather well followed by the experimental data. 
The angle of ejection a, extracted from the fits of Fig. [51 
is found to be constant for all experiments: a ~ 53° ±3°, 
which corresponds well to the corona edge angle mea- 
sured previously (Fig. Hi): 6* ~ 56° ± 3°. We have 
checked that a and 9 are robustly constant for all ex- 
periments and do not depend on any impact parameter, 
as shown in Fig.[9l No clear correlation of a with the im- 
pact energy Ec can be seen. The equality found between 
the fitted ejection angle a and the measured corona edge 
angle 9 is again in favour of the present grain ejection 
scenario. 

From eqs. (I2^)-(I2J:), we can calculate the area of the 
corona Acorit): The predicted evolution is in good agree- 
ment with experimental results and gives tA^ormax — 
1.36f//^^^ very close to the coefficient 1.38 measured ex- 
perimentally. 

Although we do not measure directly the grain ejection 
velocities, the fits of Fig. [5] give a characteristic ejection 
velocity uq- Figure llOb . shows the slow increase of uq 
with the impact velocity Uc, and its dependence on the 
sphere radius, where we can see that Uc is not the relevant 
parameter to account for the variations of uq. The kinetic 
impact energy Ec = MUc/2 is the relevant parameter, as 
shown in Fig. llOb where the kinetic energy of one ejected 
grain mu'^/2 is plotted as a function of the impact energy 
Ec in log-log axes. The best power law fit is: 
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FIG. 8: (Color online) Rescaled data of a) heights {H — 
H{0))/{H^ax-H{0)), b) top diameter _Dt/A(0) and c) bot- 
tom diameter Di,/Dt{0) as a function of rescaled times ap- 
pearing in Eqs. [3] for all of the 65 experiments. 
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FIG. 9: (Color online) Fitted ejection angle a as a function of 
the impact energy Ec for the 65 experiments (same symbols 
as in Fig. [6} . 
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FIG. 10: (Color online) Variations of ejection parameters with 
the impact parameters: a) Ejection velocity uo as a function 
of impact velocity Uc\ b) Kinetic energy of one grain mMo/2 
as a function of impact energy Ec and the best power law fit 
Eq. (|4]). Same symbols as in Fig. |6l 



FIG. 11: (Color online) a) Number of ejected grains N^j as 
a function of impact energy Ec for all experiments, and best 
power law fit given in Eq. ((Sj. b) Kinetic energy of ejecta 
_Bej as a function of impact energy Ec and best power law fit 
given in Eq. Q. Same symbols as in Fig. |6l 



m7.2/2~0.5610-X°'™', (4) 

showing a moderate effect of the impact energy on the 
kinetic energy of one ejected grain. From the ballistic 
model presented above and this experimental result of the 
scaling of the ejection velocity uq with the impact energy 
Ec, we can infer the scaling for the total duration of the 
corona ballistic dynamics T, for its maximal height Hmax 
and its maximal top diameter Dtmax with the impact 
energy Ec as: 

T = 2tH„_ = 2wosina/.goc£;0i9, (5) 
H^ax = 5iL../2cxi;0-3^ (6) 
Dtmax = ulsm2a/gxE°''\ (7) 

allowing to characterize the duration and the extension 
of deposits upon an impact. These results show that the 
corona duration T depends only slightly on the impact 
energy E^ whereas the corona expansion depends mod- 
eratly on Ec- 

In Eq. ([T]), the characteristic number of ejected grains Nej 
has been related to the maximal value of the total appar- 
ent surface of all of the ejected grains Atotmax, through 
the unknown product nw. We will estimate in section ITVl 
that nw ~ 12 mm~'^ and we use from now this value to 



calculate the number of ejected grains Ncj according to 
Eq. 11]). In Fig. [TTk is plotted as a function of the 
impact energy Ec in log-log axes. We see that whatever 
the projectile size all the experiments collapse close to 
the power law: 

iVej ~7.2105sO-"°±°-05. (8) 

We can now estimate the kinetic energy transmitted to 
all the ejected grains as Ecj — Ncjmu'^/2. Figure [Tib 
shows the evolution of Ecj as a function of Ec- The total 
kinetic energy of the ejected grains is found proportional 
to the impact energy as: 

Ecj ~ omiEc. (9) 

This unexpected result allows for t he defini tion of an ef- 
fective restitution coefficient p = y/ Ecj/Ec « 0.18. The 
small value of p confirms the well known and well used 
fact that granular beds are very efficient dissipating sys- 
tems. 



IV. DISCUSSION 

In this paper, we have shown that the number of ejected 
grains Ncj obeys a scaling law with the impact energy 
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Ec (see Eq. ([8])). In the literature, the impact energy 
has been shown to be also the relevant parameter that 
governs the crater size, in terms of its typical radius 
and depth 0, H, H, 0- Let us now relate these differ- 
ent scaling laws. The scaling laws known from Ref. Q 
for the typical radius Herat and depth Herat of craters 
arc Rcrat tx E° '^^ and Herat o<: E^-^^. The volume V of 
craters can thus be also related to the impact energy as 
V ~ HeratRlrat ^ E^'^'^ , with a numerical prefactor that 
depends on the precise shape of the crater and on the 
grain and projectile properties. From the crater profile 
measurements of Ref. [3| for glass beads, we deduce the 
precise scaling law for the crater volume: 

F ~ 3.610'^i;°-^^. (10) 

From this expression for the crater volume, we can es- 
timate the number of grains that have been ejected as 
Ncrat = V^/v where v = Airr^/S is the volume of one 
grain and $ is the solid volume fraction of the granular 
packing ($ 2± 0.6). This number of excavated grains fol- 
lows thus a power law with the impact energy with the 
exponent 0.67, that is very close to the exponent 0.70 
measured here for the ejected grains in Eq. ([5]). This def- 
initely confirms that the geometrical measurements made 
on the corona give relevant and quantitative information 
on the grain ejection dynamics. Writing that the number 
of ejected grains Nej measured here is equal to the num- 
ber of excavated grains Nerat from Q allows to bypass 
the unknow prefactor wn in Eq. ([1]) and to accurately 
calculate Nej and E^j shown in Fig. [TlJ Imposing the 
equality Nej = Nerat gives wn ~ 12 mra~'^ . 
The scaling laws uq oc E^-^^ and Nej oc obtained 
here for a large projectile over grain size ratio {R/r ~ 
25-50) are close to the results obtained for a small size 
ratio {R/r = 1) in [11]: oc E°-^'^ and WTj oc 
where uq and Nej refer to ejection velocity and number 
of ejected grains averaged over many experiments. This 
suggests that the same scaling laws hold whatever the 
ratio R/r. The constant value of the ejection angle equal 
to about 55° ±5°, is in agreement with previous values 
measured between 40° and 60° in the case of high speed 



impacts [13, H^. 

Let us recall that the energy of the ejected grains is shown 
to correspond to 3% of the impact energy (Eq. ([9])), that 
is much larger than the tiny fraction (0.3% according 
to 0) required to excavate the crater, i.e. for the grains 
to move just above the granular surface with zero veloc- 
ity. This is illustrated by the larger values of the maximal 
height of the corona Hmax (between 10 mm and 40 mm 
in Fig. 2^) than the values of crater depth (between 4 mm 
and 10 mm in [1,01) for the same material properties and 
impact conditions. This suggests to take into account the 
kinetic energy of the ejected grains when considering bal- 
ance of energies. However the essential part (^ 97%) of 
the impact energy is dissipated in the granular target 
through frictional contacts and inelastic collisions. 
V. CONCLUSION 

The dynamics of grain ejection due to a large sphere 
impacting a granular material has been experimentally 
investigated through the time evolution of the corona, 
mainly in terms of its geometry. Whereas the dimensions 
of the corona - its height, top and bottom diameters - 
change with time, the angle formed by its edge with the 
horizontal granular surface remains constant during the 
ejection process. All these geometrical properties arc well 
described by a simple ballistic model, supporting that 
grains are quasi-instantaneously ejected at time of impact 
from the same position in one direction. This direction 
appears to be constant when varying the experimental 
parameters and equal to about 55°. One may wonder 
how this angle changes with the shape of the impact- 
ing projectile. By contrast, the typical ejection velocities 
and number of ejected grains change when varying ex- 
perimental parameters and are controlled by the impact 
energy through power laws. The evaluation of the energy 
of the ejected grains allows finally for the calculation of an 
effective coefficient of restitution characterizing the com- 
plex collision process between the impacting sphere and 
the fine granular target. An important result is that this 
effective restitution coefficient is constant when varying 
the experimental parameters and equal to about 0.2. 
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